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The critical behavior of non-order parameter fields is discussed. We show that relevant features of the deconfining 
phase transition can be determined by monitoring universal properties induced by the order parameter on the physical 
excitations. Some of the behaviors we uncover are already supported by lattice results. 



1 Introduction 

Properties of phase transitions are conveniently inves- 
tigated using order parameters. However, the order 
parameter may not always be phenomenologically ac- 
cessible. Besides, in nature most fields are non-order 
parameter fields. This lead to the question whether 
it is possible to extract relevant information about, or 
even identify, the onset of a phase transition using non- 
order parameter fields. We predict a universal behav- 
ior of a non order parameter field induced by the order 
parameter field near and at the phase transition. We 
also illustrate how to extract relevant features of the 
deconfining phase transition by monitoring the critical 
properties of the physical excitations [ ^ 121 EH ED • 

To derive the general properties of a singlet field we 
couple this to the order parameter [0. In this way 
we can directly study the transfer of information, en- 
visioned in [[J, of the phase transition properties from 
the order parameter to the singlet field. The first 
studies, dealing with the information transfer from 
Polyakov loops to the hadronic states in pure Yang- 
Mills theory, are reviewed in 

Here we go beyond constant field approximation [ |2] 
and analyze two cases: time dependent order param- 
eter fields and time independent ones [|3J- Although 
time independent order parameter fields carry phys- 
ical information about the phase transition, they do 
not propagate and cannot be canonically quantized. 
An example is the Polyakov loop, which is considered 
to be the order parameter of the pure Yang-Mills the- 
ory, and by construction is a function of space only. 
Time-dependent order parameter fields are generally 
associated with physical states. These can be compos- 
ite, such as the chiral condensate, the order parame- 
ter related to chiral symmetry in QCD, or elementary, 
such as the Higgs field might be for the electroweak 
phase transition. The non-order parameter field, al- 
ways time-dependent, is a singlet in particular under 
the symmetry group whose breaking is monitored by 
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the order parameter. We assume the non-order pa- 
rameter field to have a large mass with respect to all 
the other scales in play, and hence, an associated small 
correlation length near the phase transition. A generic 
singlet field in the Yang-Mills theory is the glueball, 
which is also a physical state of the theory. 

Using general field theoretical arguments we demon- 
strate that the screening mass associated with the spa- 
tial two-point function of the singlet field is strongly 
affected by the nearby phase transition for both time- 
dependent and time-independent order parameters. 
More specifically, the screening mass has a drop near 
the phase transition. This is due to the three dimen- 
sional nature of the screening mass, which makes it 
particularly sensitive to infrared physics and it ex- 
plains the transfer of information between the order 
parameter and the singlet field. We also study the 
effects of a phase transition on the pole mass of the 
singlet field in the four dimensional theory, and find 
that only the spatial correlation lengths feel the pres- 
ence of a nearby phase transition. 

Our main conclusion in both cases is that the infor- 
mation about the phase transition, encoded in the be- 
havior of the order parameter field is transferred to, 
and obtainable from the singlet field(s) present in the 
theory. 

2 General Theory 

We consider a temperature driven phase transition and 
work in a regime close to the phase transition. In or- 
der for our results to be as universal as possible, we 
use a general renormalizable Lagrangian containing a 
field neutral under the global symmetries, and the or- 
der parameter field, as well as their interactions. The 
protagonists of our theory are two real fields, h and 
X- The field h is a scalar singlet, while \ transforms 
according to \ ~ > Z X with z S Zjy. While the gener- 
alization to Zpf is straightforward, we consider explic- 
itly the case of Z 2 , which is suitable for understanding 
the deconfining phase transition of 2 color Yang-Mills. 
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This has been heavily studied via lattice simulations 
El EJ ■ The potential is 



The full expression of the h two-point function at one- 
loop level is given by the following set of diagrams: 



V(h, X ) 
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The coefficients are real. Stability requires A > and 
<74 > 0. Assuming gi > and go < assures that the 
extremum of the potential is a minimum. This renor- 
malizable potential can be considered, for example, as 
a truncation of the one presented in [QlEli an d can be 
used to determine some of the space-time independent 
properties of the vacuum. 

We conduct our analysis using the following assump- 
tions : i) The \ field is light close to the transition, 
is massless at the transition point, hence it dominates 
the dynamics, ii) The h field is heavy, and thus we can 
neglect its quantum and Boltzman suppressed thermal 
loop corrections. 

The temperature dependence of the minimum of the 
linearized potential is qualitatively sketched in Figure 
1 of and it is discussed in [E] in more detail. 

In the following we discuss the fluctuations of h around 
its vacuum expectation value and evaluate the spatial 
correlators. 



3 3D Order Parameter 

We have chosen x = x( x ) an d decomposed the four 
dimensional h field into its Matsubara modes. After 
integrating over time the action reduces to an effec- 
tive three dimensional one. The kinetic term for h and 
its self-interaction terms receive contributions from all 
Matsubara modes. However, only the zero mode con- 
tributes to the /^-interaction, which drives the dy- 
namics of the h field close to the phase transition [ 2 , 
as the x field becomes light with respect to h. There- 
fore, we confine our discussion to the theory which 
features the fields \ and h$. For simplicity, we denote 
/io by h which is taken to be directly the fluctuation 
field around its tree level vacuum expectation value. 
The three dimensional Lagrangian reads: 



1 



1 



1 



1 



C 3 = -VWH-V X Vx+-mV + -m^ 
+ T±(x 2 ) 2 + VT 9 -±h X 2 +T^h 2 x* 



(2) 



where the coupling constants have the same mass di- 
mension as in the corresponding four dimensional the- 
ory. 
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Double lines indicate X fields and single lines stand 
for the h field. We regularize all the ultraviolet di- 
vergent contributions. The second diagram is infrared 
divergent for T — ► T c and gives the dominant one loop 
contribution to the screening mass of h. This is true 
both in the symmetric and symmetry broken phases. 
In the limit of zero external momentum 



m 2 (T) 



m 



T 
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Equation Q illustrates how the nearby phase transi- 
tion is directly felt by the non-order parameter field. 
Furthermore, it gives the general prediction that the 
screening mass of the singlet field must decrease close 
to the phase transition. If we stop the analysis here, 
at one-loop level, we predict the following critical be- 



havior for Ar 



i 2 , where m is the mass 



at a temperature close to the critical point: 



Am 2 (T) 



Am 2 (T) = - 



167TTO 

16ttM v 
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(5) 
(6) 



where M x — v2|? 



\T - T c \ y / 2 . This one-loop 



result breaks down at the transition point, due to the 
infrared singularity. But since h is not the order pa- 
rameter field, its correlation length (l/m) is not ex- 
pected to diverge at the transition point. 

We now investigate the behavior of the h screening 
mass near the phase transition by going beyond the 
one-loop approximation. 

3.1 Healing the IR behavior 

When analyzing contributions beyond one-loop order 
to the h two-point function, the number of diagrams, 
and distinct topologies one needs to consider prolif- 
erates. We select a subclass of diagrams, that heal 
the infrared divergences, while capturing the essential 
physical properties of the problem at hand. A finite 
result for the two-point function, both in the broken 
and unbroken phases will emerge, while the expression 
for the mass of h turns out to be continuous across T c . 
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In the unbroken phase it is known, that for a generic 
O(N) theory in the large N limit the following chain of 
bubble diagrams represents the leading contribution, 



These diagrams constitute a geometrical series and 
their resummation is exact [ 7 . We choose the same 
set of diagrams in the unbroken phase. We note that 
a large N approximation might not be the best choice 
for investigating the details of the phase transition, 
however we show that it well reproduces the behav- 
ior of the screening mass of hadronic degrees of free- 
dom near the deconfinement phase transition as given 
by lattice simulations. After summation the following 
expression for m follows: 



m 2 (T) 



T 



9l 



167TTO V + XT 



(7) 



which is finite at T c , where it yields: 



m 2 (T c ) = to 2 



9l 



(8) 



This is the main result of [E], predicting that close 
to the phase transition the singlet state must have 
a decreasing mass parameter, associated with spatial 
correlations. More specifically, the drop at the phase 
transition point is given by the ratio of the square of 
the coupling constant governing the interaction of the 
singlet state with the order parameter (gi) to the or- 
der parameter field self-interaction coupling constant 
A. In this way, via the drop of the singlet field at 
the phase transition, one can derive a great deal of in- 
formation about the phase transition, and about the 
order parameter itself. 

The analysis is much more complicated in the broken 
phase. Here \ develops an average value, v, which 
induces one also for h [|3J. So on general grounds, the 
fields \ and h mix in this phase. The mixing angle is 
proportional to giv/m 2 , and therefore, this mixing 
can be neglected within the present approximations. 
Like for the symmetric phase, we consider only the 
effects due to the x loops for the h propagator. Due 
to symmetry breaking we now need to include also the 
trilinear \ coupling 



A 
3! 



(9) 



which is expected to affect the analysis. Curing the 
infrared divergence is thus now more involved due to 
symmetry breaking. In the case of the large N limit 



of O(N) symmetry, one can show, that diagrams with 
trilinear vertices are again suppressed relative to the 
simple bubble diagrams. But now we go beyond the 
large N limit by computing a new set of diagrams, 
which can be evaluated exactly, and thus capture rele- 
vant corrections due to symmetry breaking, neglected 
in the large N limit. The new chain of diagrams we 
compute has terms of the form: 



-00 • • 00 • • <D- 



(10) 



This class of diagrams has knowledge about the onset 
of symmetry breaking via the presence of the trilinear 
vertices, and, in that respect, is the simplest exten- 
sion of the chain of simple bubble diagrams. Another 
amusing property of 1)10(1 is, that the sum can be per- 
formed exactly, yielding again an infrared finite result: 



m 2 (T) 



24irM Y + XT 



X 2 T 2 + 2AnXTM y 



384ir 2 M 2 



(11) 



We see that m(T c ) from the broken side of the tran- 
sition equals exactly the one from the unbroken side 
of the transition, even when departing from the large 
N limit. The mass squared of h is a continuous func- 
tion through the phase transition, and the associated 
correlation length remains finite. This result does not 
hold order by order in the loop expansion, but only 
when the infinite sum of the diagrams is performed. 

In order to disentangle relevant properties of the phase 
transition, we construct slope parameters for the sin- 
glet field: 



1 dm 2 (T) 
hm — — — — 

n^ T ± Am 2 (T) 



dT 



(12) 



with Am 2 (T c ) = g\jX . The functional form of T> + 
and T>~ is the same, provided that the same class of 
diagrams is resummed on both sides of the transition. 
In such case: 



XT 



16 7r ,. dm Y 
- hm — - 

X T c t->t- dT 



(13) 



in the symmetric phase, and T> + in the broken phase is 
obtained by replacing m x with M x = >/2\m x \. While 
the mass of h remains finite at T c , its slope encodes 
the critical behavior of the theory. For example, if 
m 2 vanishes as (T c — T) v close to the phase transition 
(with the correlation length £ cx \T — T c \~ v / 2 ), then 
scale with exponent [y/2 — 1). A difference in the 
functional form of the slopes may emerge, when on the 
two sides of the transition different class of diagrams 
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Figure 1. Left panel: Screening mass of the singlet field 
close to the phase transition as function of the tempera- 
ture. Right panel: Lattice data; figure from 

are resummed. Using the wider class of diagrams in 
the broken phase considered above while retaining just 
the simple bubble sum in the unbroken we determine: 



(14) 



This is due to the onset of spontaneous symmetry 
breaking, i.e. in the broken phase the class of re- 
summed diagrams contains trilinear type of interac- 
tions. We identify thus a less singular behavior with 
respect to the simple sum of bubbles. More specif- 
ically, the scaling exponent for T> + is now {y — 1). 
While the explicit relations between the scaling expo- 
nents and the slopes are only valid within the given 
summation scheme, these quantities are, nevertheless, 
a measure of the critical behavior near the phase tran- 
sition. Only experimental results will be able to select 
which class best describes the data. 

In figure^we schematically present the behavior of the 
h screening mass as a function of temperature, in units 
of the critical temperature, for m 2 oc (T c — T) . The 
left panel illustrates the rapid decrease of the singlet 
field screening mass in the critical region. The right 
panel is the lattice data of [ 8 obtained for the glueball 
screening mass in the three color Yang-Mills theory. 
The resemblance between our results and the lattice 
results is intriguing, but more lattice data is needed 
in order to distinguish between the above possibilities, 
and quantitatively determine the size of the drop. 

4 4D Order Parameter 

Here we present the effects of the four dimensional 
X — x( x j t) order parameter field on the properties of 
a singlet h — h(x, t) field. The order parameter field 
is now a physical field, and as such, it can propagate 
in time, and can be canonically quantized. The four- 
dimensional Lagrangian of the renormalizable theory 
we use to define our Feynman rules is: 



1 



1 



-h z 



?x 2 



4!* 



2± h x 2 - XL h 2 x 2 - ^h 3 ~ ^h\l5) 
2 4 3! 4! v ' 



All coupling constants are real and A, #4 > 0. The rele- 
vant and marginal couplings, in terms of dimensionless 
ones, are g\ = g\m , 53 = g^ra and 52 = 52 ■ The 
X field is subject to Zi symmetry. In four dimensions 
a Bose-Einstein distribution emerges for \, and thus 
thermal fluctuations become important. These have 
the tendency to restore symmetry, destroying the pos- 
sibility for the formation and existence of any physical 
condensate at high temperatures. 

The diagrams contributing to the h two-point func- 
tion at one- loop order are shown in Q. Even though 
these are the same as in the three dimensional theory, 
here they are computed taking into account the time 
dependence of x- We can show, that finite tempera- 
ture corrections from all diagrams that involve an h 
loop are Boltzman suppressed. Thus the relevant dia- 
grams for the singlet h field are those on the first line 
in JSJ • Unlike the three dimensional case where ultra- 
violet divergent tadpoles have been removed by renor- 
malization, in the four dimensions these diagrams also 
provide temperature contributions. The diagrams are 
evaluated using standard techniques of the imaginary 
time formalism The tadpoles are real and provide 
temperature dependence to the h mass, not present 
in the three dimensional theory. The eye diagram 
has a real and an imaginary part contributing to 
the two-point function [E 2 — m 2 — II(J5)] 1 . Here 
m is the tree level mass of h and at one-loop order 
n = Utadpole + ^eye ■ For zero external momentum 



TL eye (E) =- — 



d 3 k 1 



(2tt) 3 u E 2 - Alu 2 



(1+ 2/( W )),(r6) 



with u) = \J\t 2 + rn x 2 , and / Bose distribution func- 
tion. The real part is given by the principal value 
of (fTB)l . and represents a shift in the mass squared of 
the h. For h at rest, the real part by its definition 
determines the pole mass M, i.e. the pole of the full 
two-point function 



M 2 



n. 



tadpole 



Re U eye (E — M) = . 



(17) 



Self-consistent solution of the above equation shows 
that the large tree level mass dominates, and loop cor- 
rections are negligibly small in the temperature range 
of interest, near the phase transition. Thus M ~ m 
acts as an infrared cutoff, guaranteeing the absence of 
infrared divergence in the pole mass. All the one-loop 
contributions to the pole mass sum to 



M 



1 + (92 - hh - 2.9i 1 



T 2 



24 m 2 



(18) 
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It is important to distinguish between the pole mass 
and the screening mass. The screening mass, m s , is 
defined by the location of the pole in the static prop- 
agator for complex momentum p — im s , 



P 



T1(E = 0,p) = 0. 



(19) 



In the small momentum limit this leads to the follow- 
ing definition 



ttT + lim Tl(E 

p^0 



0,p). 



(20) 



In three dimensions the pole and screening masses are 
the same, since the static propagator is the full propa- 
gator. We have shown, that when looking at the pole 
mass the IR problem of the eye diagram is regulated 
by the heavy h mass. When analyzing the screening 
mass, however, the one-loop IR divergence present in 
the three dimensional case is recovered. This can be 
seen by setting E — in expression (|16l) corresponds 
to reducing the 4D theory to a 3D one. To display the 
relevant infrared contribution we take the high tem- 
perature expansion which yields: 



v 



T 



(p 2 + to 2 ) 2 



32-^771} 



(21) 



Thus, for the static limit the phase transition region is 
dominated by the same type of infrared divergence we 
encountered in the time-independent order parameter 
field case. Note, however, that the numerical constant 
in eq. (|21() differs from that in the second term of eq. 

The reason for this is that the reduction was 
done only for the modes of h, and thus in eq. (|21|1 . 
contrary to (0J, all Matsubara modes of x contribute. 
By combining all the tadpole and eye diagrams we find 
for the screening mass at one-loop order 



, 9l T T 2 
"(Sift "ft) 



32-7T m-i 



(22) 



showing clearly the eye contribution as the infrared 
dominant one. Note also, that above we have tree 
level coefficients, and a complete investigation would 
require renormalization group analysis. 

5 Deconfinement and Conclusions 

When analyzing many physical situations, it is com- 
mon practice to isolate the order parameter field, and 
more generally, the light degrees of freedom, since 
these are expected to be the relevant states at low 
energies. While this procedure certainly is reasonable, 
in nature most of the physical fields are neither order 



parameter fields, nor light at all. In order to extract 
information from these heavy states, we needed first to 
determine new and universal features associated with 
them. We have used a general strategy proposed first 
in [HE]; according to which we couple the light de- 
grees/order parameter fields to the heavy fields in the 
most general way, and then truncate the theory by re- 
taining all the relevant and marginal operators in the 
Lagrangian. In doing so, the theory is fully renormal- 
izable while capturing the relevant contributions. The 
operator set is further constrained by imposing the rel- 
evant symmetries of the problem at hand. In order for 
our procedure to work, we also assume, that the other 
physical states of the theory have masses larger than 
our non-order parameter field. In this way we can, 
formally, integrate these states out, and their effects 
are absorbed in the modified couplings of our effective 
Lagrangian. 

For both the time independent and time dependent 
order parameter field we have shown [ [3], that the 
spatial correlators of the non-order parameter field arc 
infrared dominated, and hence can be used to deter- 
mine the properties of the phase transition. We have 
determined the general behavior of the screening mass 
of a generic singlet field, and have shown how to ex- 
tract all the relevant information from such a quantity. 
We have further demonstrated, that the pole mass of 
any non-order parameter physical field is not infrared 
dominated. 

Our results can be immediately applied to any generic 
phase transition. We have used as relevant example, 
for the time independent order parameter field case, 
the deconfining transition of Yang-Mills theories. In [ 
^ a model containing glueball H and Polyakov loop 
£(x) was proposed, and is supported by other investi- 
gations The present results can be understood 
as the higher loop corrections to the glueball model in 
and can be immediately applied to the two color 
Yang-Mills phase transition, once we associate H with 
h, and I with \- For example one can take the follow- 
ing relation between H and the glueball field h: 



H=(H)(l + 



(23) 



Here (H) = A 4 /e is the vacuum expectation value of 
the glueball field below the critical temperature and c 
is a positive dimensionless constant fixed by the mass 
of the glueball. For the £ field we have \ — with 
k a mass dimension two constant, which at high tem- 
perature is proportional to T 2 . Specifically, our renor- 
malizable theory is a truncated (up to fourth order in 
the fields) version of the full glueball theory. 

We confronted already our theoretical results with lat- 
tice computations of the glueball screening mass be- 
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havior close to the phase transition studied in [ |H] 12. A. Mocsy, F. Sannino and K. Tuominen, 
for three colors. Our analysis not only is in agree- Confinement versus chiral symmetry, 

ment with the numerical analysis, but also allows us arXiv:hep-ph/0308135 
to provide a better understanding of the physics in 
play. Further lattice simulations are able to deter- 
mine the coupling strength of any glueball state to 
the Polyakov loop by following the temperature de- 
pendence of screening masses of such states. 

Our analysis suggests, that monitoring a number of 
spatial correlators, or more specifically their deriva- 
tives, is an efficient and sufficient way to experimen- 
tally uncover the chiral/deconfining phase transition 
and its features. 

The findings here provide a general computational 
strategy useful to deduce new quantitative information 
about any phase transition. Using this strategy we 
were also able to explain how deconfincmcnt is linked 
to chiral symmetry restoration [I12|. 
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